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ABSTRACT 

Due to a new global analysis method, it is now possible to measure the internal composition of pulsating white 
dwarf stars, even with relatively simple theoretical models. The precise internal mixture of carbon and oxygen 
is the largest single source of uncertainty in ages derived from white dwarf cosmochronometry, and contains 
information about the rate of the astrophysically important, but experimentally uncertain, 12 C(a,7) 16 nuclear 
reaction. Recent determinations of the internal composition and structure of two helium-atmosphere variable 
(DBV) white dwarf stars, GD 358 and CBS 114, initially led to conflicting implied rates for the 12 C(a,7) 16 
reaction. If both stars were formed through single-star evolution, then the initial analyses of their pulsation 
frequencies must have differed in some systematic way. I present improved fits to the two sets of pulsation 
data, resolving the tension between the initial results and leading to a value for the 12 C(a, 7) 16 reaction rate that 
is consistent with recent laboratory measurements. 
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1. INTRODUCTION 

As a white dwarf star cools over time, from a hot planetary 
nebula nucleus down to the coolest temperatures possible given 
the finite age of the universe, there are three narrow ranges of 
surface temperature where they may become pulsationally ex- 
cited. The particular interval where a given white dwarf will 
pulsate is determined by its spectral type. The PG 1 159 stars or 
variable DO white dwarfs are the hottest class, and require de- 
tailed evolutionary models for an accurate pulsational analysis. 
The cooler helium-atmosphere (DB) and hydrogen-atmosphere 
(DA) variables are both easier to model, since relatively simple 
polytropic models may be used to approximate their evolution 
quite accurately (Wood 1990). The DB variables (DBVs) are 
structurally the simplest, since they have a single surface layer 
of helium on top of a presumed carbon/oxygen (C/O) core, 
while DA variables (DAVs) have an additional surface layer of 
hydrogen above this helium mantle. 

Although our theoretical models of variable white dwarfs 
are relatively simple, they have been able to match the ob- 
served pulsation periods of these stars with a degree of accu- 
racy that — with the exception of the Sun — is unsurpassed in 
the field of stellar seismology. In part, this achievement has 
been made possible by the simplifying physical circumstances: 
the absence of nuclear fusion, the extremely high surface grav- 
ity (logg ~ 8), slow rotation (P rot ~ 1 d), and negligible mag- 
netic fields. However, the most recent improvements in the 
agreement between models and observations have been driven 
more by global explorations of the defining parameters of exist- 
ing models, rather than by fundamental changes to the models 
themselves. 

By combining improved models with objective global opti- 
mization methods, we may eventually obtain physically com- 
plete and fully self-consistent theoretical models of white dwarf 
stars. Until then, we can rely on the relative quality of the match 
between our models and observations as a proxy for the abso- 
lute optimization that would be possible if our models were 
physically complete descriptions of real white dwarfs. This 
may lead to systematic errors in the values of some of our pa- 



rameters, depending on the nature of the incompleteness, but it 
will allow us to proceed and to develop useful algorithms for 
diagnosing the limiting assumptions in our models. 

Recently, Metcalfe & Charbonneau (2003) have given a de- 
tailed description of a method developed to match simple mod- 
els of DBV white dwarfs to the available observations. The 
method itself is perfectly general, and may easily be extended 
to other types of pulsating stars. In the context of DBV 
white dwarfs it has been used, among other things, to infer 
the internal mixture of carbon and oxygen — the most impor- 
tant source of uncertainty in age estimates from white dwarf 
cosmochronometry — and, by extension, the rate of the key 
12 C(a,7) 16 reaction. The initial applications of this method to 
two white dwarfs, first to the brightest DBV star, GD 358 (Met- 
calfe, Salaris, & Winget 2002), and then to the faintest mem- 
ber of the class, CBS 1 14 (Handler, Metcalfe, & Wood 2002), 
yielded apparently conflicting results. Here, I will demonstrate 
that the apparent conflict arose from systematic differences in 
how the two data sets were analyzed. When the optimum mod- 
els are determined using identical criteria, both stars yield in- 
ternal compositions that imply rates for the 12 C(a,7) 16 reac- 
tion that are consistent with one another, and which both agree 
with recent extrapolations from high-energy laboratory mea- 
surements. 

2. INITIAL RESULTS 

The global optimization method for DBV white dwarf mod- 
els, based on the publicly available genetic algorithm PIKAIA 
(Charbonneau 1995), was developed by Metcalfe, Nather, & 
Winget (2000) and later extended by Metcalfe, Winget, & Char- 
bonneau (2001, hereafter MWC) to allow five adjustable pa- 
rameters. The method attempts to minimize the differences be- 
tween the observed and calculated periods and period spacings 
for models with effective temperatures (T e s) between 20,000 
and 30,000 K, total stellar masses (M*) between 0.45 and 0.95 
M©, helium layer masses with -log(MH e /A^*) between 2.0 and 
~7.0, and an internal C/O profile with a constant oxygen mass 
fraction (Xq) out to some fractional mass (q) where it then de- 
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creases linearly in mass to zero oxygen at 0.95 m/M*. The 
broad range of the search is limited only by observational con- 
straints and by the physics of the models. 

The observational data for GD 358 came from an extremely 
successful multi-site campaign in 1990 by the Whole Earth 
Telescope collaboration (WET; Nather et al. 1990). The results 
and analysis of these observations were reported by Winget et 
al. (1994) and Bradley & Winget (1994, hereafter BW) who 
conclusively identified a series of 1 1 consecutive radial over- 
tones (A:=8-18) of non-radial g-mode pulsations with the same 
spherical degree (1=1). Using this rich data set, BW attempted 
to match the periods themselves (Pj) as well as the spacings 
between consecutive modes (AP = P^+i —Pk). When confined 
to search the same range of parameters considered in this initial 
study, MWC found the same optimal set of model parameters as 
BW. But the global search revealed a significantly better model 
outside the range of the initial search, with optimal parame- 
ters r eff = 22,600 K, M* = 0.650 M , log(M He /M«) = -2.1 A, 
Xo = 0.84, q = 0.49 m/M*, and root-mean-square (rms) residu- 
als ap = l .28 s, and cta/> = 1 -42 s. 

The central oxygen mass fraction in white dwarf models is 
primarily determined by the rate of the 12 C(a,7) 16 reaction. 
The extent and efficiency of internal mixing during the red gi- 
ant phase is a secondary factor (see §4). This allowed Metcalfe, 
Salaris, & Winget (2002) to adjust the n C(a, 7) lfi O rate in evo- 
lutionary models of the internal chemical profiles until they pro- 
duced the optimal value of Xo for GD 358. The implied reac- 
tion rate was S300 = 370 ± 40 keV barns, significantly higher 
than most extrapolations from high-energy laboratory measure- 
ments. The authors noted that the application of this method to 
additional DBV white dwarfs could provide independent mea- 
surements of this important nuclear reaction rate. 

Although CBS 114 is the faintest known DBV, single-site 
observations from SAAO in 2001 (Handler, Metcalfe, & Wood 
2002), combined with a reanalysis of the discovery data of 
Winget & Claver (1988), revealed 7 stable pulsation modes 
with periods between about 400-650 seconds. The mean period 
spacing was consistent with non-radial g-mode pulsations of 
spherical degree i=\, and a comparison with models identified 
the radial overtones as fe=8,9 and 11-15, with the k= 10 mode un- 
detected. The global analysis of these pulsation data had to be 
restricted to comparing the observed and calculated periods — 
neglecting the period spacings because of the undetected mode 
that interrupts the sequence of otherwise consecutive radial 
overtones. The optimal values for the five model parameters 
were r eff = 21,000 K, M* = 0.730 M Q , log(M He /M*) = -6.66, 
Xo = 0.61, q = 0.51 m/M*, and the rms period residuals were 
up = 0.43 s. This fit matched the observed periods so closely 
that it also reproduced the observed deviations from the mean 
period spacing (an alternate period spacing criterion) almost ex- 
actly (see Handler, Metcalfe, & Wood 2002, Fig. 9c). 

The higher stellar mass for CBS 1 14 relative to GD 358 
should naturally lead to a slightly lower central oxygen mass 
fraction, since the 12 C(a,7) 16 reaction is slightly less effi- 
cient at higher densities. However, the internal chemical pro- 
file models with the appropriate mass required a rate for the 
12 C(a,7) 16 reaction of only S300 ~ 180 keV barns to repro- 
duce the optimal central oxygen mass fraction. This is very 
close to the rate derived from recent laboratory measurements 
(S300 = 165 ± 50 keV barns; Kunz et al. 2001). 

The conflicting results for GD 358 and CBS 1 14 led Handler, 
Metcalfe, & Wood (2002) to speculate that either some source 



of systematic uncertainty in the models was affecting the anal- 
ysis of the two stars in different ways, or that they may have 
different evolutionary origins. The easiest way to test the sec- 
ond idea would be to apply the method to an additional DBV 
star, and see if the internal composition agreed with either of the 
other two. Although the WET recently conducted a multi-site 
campaign on the DBV star PG 1456+103, the analysis of those 
data are still in progress. Thus, I have attempted to follow up on 
the first possibility, and identify possible sources of systematic 
error. 

One obvious systematic problem with both of the fits, as 
noted in each of the initial analyses, was the disagreement be- 
tween the optimal masses and temperatures and the values of 
these parameters from the spectroscopic analysis of Beauchamp 
et al. (1999). Both global fits produced masses that were sys- 
tematically high, and temperatures that were systematically 
low. Some discrepancy was expected, since the two methods as- 
sumed slightly different values for the mixing-length/pressure 
scale height ratio for atmospheric convection. In addition, the 
simple polytropic models used for the global fits included static 
envelopes containing 5% of the total stellar mass. Metcalfe, 
Salaris, & Winget (2002) demonstrated that by reducing the 
fractional mass of these static envelopes to only 2%, the mass 
and temperature of the optimal model for GD 358 were in closer 
agreement with the spectroscopic values. In any case, neither 
of these two problems should be responsible for the differences 
between the results for GD 358 and CBS 1 14, since both prob- 
lems were common to the two analyses. 

Upon reflection, there were two things that did differ between 
the analyses of the two stars: (1) the fit to GD 358 used the pe- 
riods and the period spacings to judge the quality of the match, 
and (2) there were a larger number of higher radial overtone 
modes in the case of GD 358. To evaluate whether or not these 
subtleties could be responsible for the differences between the 
derived values of the central oxygen mass fraction, we can look 
at a small grid of models with various values of Xq and q keep- 
ing the mass, temperature, and helium layer mass fixed at their 
optimal values. The effect of these two differences on the de- 
termination of the central oxygen mass fraction for GD 358 are 
illustrated in Figure 1 . Each panel shows the optimal combi- 
nation of Xo and q for a given fitting criterion in black, along 
with the ranges of model parameters leading to rms residuals 
that differ from those of the optimal combination by 1, 3, and 
10 times the observational uncertainty (cr bs ~ 0.03 s), shown 
in progressively lighter shades of grey. The difference between 
the distributions in the top two panels demonstrates that when 
the period spacings are used in addition to the periods them- 
selves for fitting, the effect is to bias the central oxygen mass 
fraction toward higher values and to increase the uncertainty in 
the optimal value of Xo- Comparing the top and bottom pan- 
els, it is clear that attempts to match only those radial overtones 
observed in CBS 114 may lead to systematically lower values 
for Xq. In principle, these two results would lead us to expect 
that the initial determination of Xo for GD 358 may have over- 
estimated the actual value due to the use of period spacings for 
fitting, and the value for CBS 1 14 may have been an underesti- 
mate since there were fewer modes of higher radial overtone. In 
practice, this exercise cannot provide quantitative information 
about the expected shifts, since the three model parameters that 
were fixed for Figure 1 will undoubtedly shift slightly in a full 
re-optimization. 
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FIG. 1 . — The optimal combinations of Xq and q for GD 358 (black) and 
the distribution of parameter values that yield rms residuals within 1, 3, and 
1 times the observational uncertainty (progressively lighter shades of grey) 
when the mass, temperature and helium layer mass are fixed at their optimal 
values from MWC and when the fitting criterion is based on matching (a) the 
pulsation periods only, (b) the period spacings only, and (c) the periods of only 
those pulsation modes observed in CBS 1 14. 

3. NEW MODEL-FITTING 

Motivated by the qualitative evidence from Figure 1, I per- 
formed several new fits using the global model-fitting proce- 
dure of MWC. Since I also wanted to reconcile the discrep- 
ancies between the masses and temperatures inferred from as- 
teroseismology and spectroscopy, the new fits used the so- 
called "ML2/a=1.25" prescription for convection adopted for 
the spectroscopic fits by Beauchamp et al. (1999) rather than the 
"ML3" prescription used by MWC. Both prescriptions follow 
the same mixing-length theory of Bohm & Cassinelli (1971), 
but assign different values (1.25 and 2.0 respectively) to the 
mixing-length/pressure scale height ratio. Each of the new fits 
use only the periods (and not the period spacings) to determine 
the optimal model, which also makes them directly comparable 
to the alternative model fit of Fontaine & Brassard (2002). 

The values of the globally optimal model parameters for the 
new fits are listed in Table 1, along with the rms differences be- 
tween the observed and calculated periods (erp) 1 . Note that 
the period residuals for the "GD 358 [11]" fit are consider- 
ably smaller than in the fit of MWC (which had a P = 1.28 s). 
The period residuals for CBS 1 14 also improved slightly, com- 
pared to the fit of Handler, Metcalfe, & Wood (2002, which 
had (7 p = 0.43 s). Small grids of models with various combina- 
tions of Xq and q, keeping the other parameters fixed, revealed 
that the formal la statistical uncertainties on the central oxygen 
mass fraction for GD 358 and CBS 1 14 are near ±0.01 in both 
cases, but I adopt an uncertainty twice this size. 

Using the mass and central oxygen mass fraction to de- 
rive a measurement of the 12 C(a,7) 16 reaction rate is model- 
dependent. To facilitate comparison with previous results, I re- 
lied on the internal chemical profiles of Metcalfe, Salaris, & 

' An additional fit for GD 358, using ML2/a=1.25 convection and both P and 
o~p= 1.32 and <ta/>=1-50. This suggests that the lower value for Xq is primarily 
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FIG. 2. — The ±lcr range of theoretical internal chemical profiles for a 0.65 
Mq white dwarf model using the NACRE rate for the 12 C(a,7) l6 reaction 
(Angulo et al. 1999, shaded area), along with profiles scaled to the optimal cen- 
tral oxygen mass fraction derived for GD 358 (upper hashed area) and CBS 114 
(lower hashed area) with the corresponding values for the 12 C(o,7) 16 rate. 

Winget (2002) and M. Salaris (private communication). In- 
terpolating between these models, the 12 C(a,7) 16 reaction 
rate that was required to match the optimal values of the mass 
and central oxygen mass fraction for GD 358 and CBS 114 
was S300 = 195 ±15 keV barns and S300 = 190 ± 15 keV barns 
respectively. Both measurements are in excellent agreement 
with the rate suggested by the NACRE collaboration (S300 = 
200 ± 80; Angulo et al. 1999). This remarkable agreement is il- 
lustrated in Figure 2, which shows the ranges of internal chem- 
ical profiles corresponding to the NACRE rate for a 0.65 Mq 
white dwarf model (from Metcalfe, Salaris, & Winget 2002) 
along with profiles scaled to match the optimal values of the 
central oxygen mass fraction for GD 358 and CBS 114. Note 
that the I2 C(a, 7) I6 rates derived from the two stars are statis- 
tically indistinguishable, but the chemical profiles themselves 
do not overlap at the center because the models have differ- 
ent masses. Both measurements are also consistent with the 
more recent laboratory value derived by Kunz et al. (2001, 
S 30 o= 165 ±50). 



4. DISCUSSION 

A re-analysis of the pulsation data for GD 358 and CBS 1 14, 
with ML2/a=1.25 convection and using only the periods to de- 
termine the optimal models, yielded improved fits leading to 
measurements of the 12 C(a,7) 16 reaction rate that are both 
consistent with each other and with laboratory measurements. 
Since a larger number of pulsation modes are present in GD 358 
to constrain the fit, it should be considered the more reliable 
measurement. The influence of having fewer modes and lower 
radial overtones on the derived central oxygen mass fraction is 
ambiguous: for a given mass, temperature, and helium layer 
mass, it may lead to an underestimate — but when all parame- 



AP led to the optimal parameters (22,900 K, 0.645 Mq, 
due to the change of convective prescription. 



-2.69, 0.63, 0.46), yielding 
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Table 1 

New model-fits for two DBV stars 



Parameter GD 358 [1 1] GD 358 [7] CBS 1 14 

T eff (K) 22,900 22,100 20,300 

M* (AT©) 0.660 0.665 0.750 

log[Af He /Af»] . . . -2.00 -2.79 -6.82 

X 0.67 0.76 0.63 

q(m/M*) 0.48 0.47 0.51 

a P (s) 1.05 0.52 0.41 



ters are allowed to re-adjust (as in the "GD 358 [7]" fit in Table 
1) it may lead to an overestimate. This underscores the impor- 
tance of a global analysis when trying to determine the optimal 
model. 

The derived central oxygen mass fractions for these two 
DBVs are now perfectly consistent with each other — with a 
slightly lower value for the more massive white dwarf, as ex- 
pected. Mapping the internal chemical composition into a rate 
for the 12 C(a,7) le O reaction is model-dependent, as discussed 
recently by Straniero et al. (2003), because the adopted mix- 
ing scheme and numerical treatment of breathing pulses are 
both important to the outcome. Unless the detailed shape of 
white dwarf internal chemical profiles can be measured empiri- 
cally through asteroseismology, or until some independent con- 
straints on internal mixing are possible, this systematic uncer- 
tainty will remain. 

Switching to the same convective prescription as was used in 
the spectroscopic analysis of Beauchamp et al. (1999) did not 
resolve the systematic errors in the derived mass and tempera- 
ture of either object. This outcome appears to be related to the 
simplified treatment of the surface helium layer and the use of 
static envelopes. 

Models that include time-dependent diffusion in the en- 
velopes, such as those of Dehner & Kawaler (1995), reveal that 
the surface helium layers of DBV stars may contain a double- 
layered structure — with a relatively thin pure He layer overly- 
ing a thicker mixed He/C/O layer above the C/O core. The 
most extensive observational test for this structure has, so far, 
come from Fontaine & Brassard (2002) who calculated a grid 
of models with masses and temperatures near the spectroscopic 
values for GD 358 and compared the observed and calculated 
periods. By including these double-layered envelopes — but ne- 
glecting the theoretically expected composition and structure in 
the core — the authors were able to match the pulsation periods 
of GD 358 with the same level of precision as the fit of MWC, 
which used single-layered envelopes, ML3 convection, and a 
physically motivated C/O core. 

An analytical explanation of how these two models — one 
with extra structure in the core, and the other with extra struc- 
ture in the envelope — can produce fits of comparable quality 



was recently provided by Montgomery (2003). Composition 
transition zones at certain locations in the core and envelope 
can produce deviations from uniform period spacing that are 
difficult to distinguish from each other. This leads to a potential 
ambiguity between, for example, the expected C/O transition 
near 0.5 m/M* and the outer He transition near 10~ 6 m/M*. As 
a consequence, it is crucial that we determine which feature has 
the largest imprint on the observed periods. 

The new fit for GD 358 in Table 1 is more directly compa- 
rable to Fontaine & Brassard's fit, since it also uses only the 
observed periods to determine the optimal model. Even after 
accounting for the additional free parameter, the new fit has 
significantly lower residuals than the pure C double-layered en- 
velope fit, suggesting that the internal C/O profile may be the 
more important of the two internal structures. However, the 
continued presence of a local minimum with log(M He /A^*) ~ 
-6.0 in fits using single-layered models suggests that double- 
layered envelopes may improve the fit even further. 

Fontaine & Brassard have also concluded that improved 
fits to GD 358 may be possible using models that include 
both double-layered envelopes and a variable core composi- 
tion. However, their calculations were limited to model grids 
using several uniform C/O mixtures, so they did not have access 
to any information contained in the residuals that might allow 
a determination of the C/O profile. In the initial results from 
an ongoing study (Metcalfe et al., in preparation) we find that 
double-layered envelope models with an adjustable C/O profile 
yield significantly improved fits relative to models with a pure 
C core. 

Clearly, the way forward is a thorough exploration of hy- 
brid models with self-consistent double-layered envelopes and 
physically-motivated internal C/O profiles. With some luck, 
this may eventually reduce the rms residuals down to the level 
of the observational noise, and provide the most accurate value 
possible for the 12 C(a, 7) ie O reaction rate. 

I would like to thank Ed Nather, Don Winget, Paul Charbon- 
neau, Maurizio Salaris, and Gerald Handler for fruitful collab- 
orations during the various stages of this project. 
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